Steady state entanglement between hybrid light-matter qubits 

Dimitris G. Angelakis^'^P Stcfano Mancini'^[!j and Sougato Bosc tH 

^ Centre for Quantum Computation, Department of Applied Mathematics and Theoretical Physics, 

University of Cambridge, Wilberforce Road, CB3 OWA, UK 

'^Science Department, Technical University of Crete, Chania, Creece 73100 

^ Dipartimento di Fisica, Universitd di Camerino, 1-62032 Camerino, Italy and 

"^Department of Physics and Astronomy, University College London, Cower St., London WCIE 6BT, UK 

(Dated: May 9, 2008) 

We study the case of two polaritonic qubits localized in two separate cavities coupled by a 
fiber/additional cavity. We show that surprisingly enough, even a coherent classical pump in the 
intermediate cavity/fiber can lead to the creation of entanglement between the two ends in the 
OO ■ steady state. The stationary nature of this entanglement and its survival under dissipation opens 

^«— ^ ' possibilities for its production under realistic laboratory conditions. To facilitate the verification of 

^>—^ , the entanglement in an experiment we also construct the relevant entanglement witness measurable 

^N . by accessing only a few local variables of each polaritonic qubit. 

C^ ^ PACS numbers: 03.67.Mn, 42.50.Ct, 03.65.Yz 
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I. INTRODUCTION 



Recently, there has been a growing interest in exploiting a certain class of coupled hybrid light-matter systems, 

namely coupled cavity polaritonic systems, for various purposes such as for realizing schemes for quantum computation 

[l|, , for communication Q and for simulations of quantum many-body systems [j, y, 0, 0, [^, [^, [lO, [ll| • These cavity- 

^ ' atom polaritonic excitations are different from propagating polaritonic excitations in atomic gases and exciton-photon 

O^. polaritons in solid state systems ^^- This area is also distinct from those using hybrid light-matter systems in 

quantum computing where only the matter system (such as an atom or an electron) acts as the qubit. In the latter 

fvj i case the qubits are atoms and light is used exclusively as a connection bus between them [ij, [ij, [la, Mj [l3j H, ^M ■ 

^ ' Promising schemes to produce steady state entanglement between atoms in distinct cavities have also been proposed 

C^ I [3 ■ III these ground states of atoms have been used in order to circumvent decoherence due to spontaneous emission. 

C^ . In addition to auxiliary atomic levels, external driving fields as well as an unidirectional coupling between cavities 

Q^ ' are required. In polaritonic coupled cavity systems on the other hand, the localized mixed light-matter excitations, 

or polaritons, allow for the identification of qubits that possess the easy manipulability and measurability of atomic 

qubits, while also being able to naturally interact whereas separated by distances over which photons can be exchanged 

between them. Motivated by the rapid experimental progress in Cavity Quantum Electrodynamics and the ability 

^^ ' to couple distinct cavities in a variety of systems [20, lU, [23, [2^, [2J| , the realization of a system that could produce 

verifiable, steady state entanglement between two polaritonic qubits in currently realistic laboratory conditions would 

be extremely interesting. In that case the decoherence emerging from the photonic losses due to the mixed nature of 

the polaritons, in addition to that from atomic spontaneous emission, will need to be controlled. Therefore, apriori 

one may not expect a completely stationary entanglement of two polaritons unless the unavoidable loss of coherence 

C^ ' due to both channels can somehow be "re-injected" into the system. 

Here we show that even under strong dissipation in both the atomic and photonic parts, it is still possible to 
detcrministically entangle two such polaritonic qubits. More precisely, we study the case of two polaritonic qubits 
coupled by a fiber/additional cavity and show that surprisingly enough, even a coherent classical pump can lead to 
the creation of entanglement between them in the steady state. The stationary nature of this entanglement should 
make easier its experimental verification. To this end we also provide a relevant operator (an "entanglement witness" 
[25| ) measurable by only measuring local variables of each polariton. 
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II. THE MODEL 

The Haniiltonian describing an array of N identical atom-cavity systems is the sum of the free hght and dopant 
parts and the internal photon and dopant couplings 

N N 



if/'-- = ^d^4afc+L^o^|e)fc(e|, (1) 

fe=i fc=i 

N 

H"'' = 9Y,{ai\g),{e\+a,\e)k{g\). (2) 

k=l 

Here a^, a^ are the photonic field operators localized in the fc-th system and |e)fc, \g)h are the excited and ground state 
of the dopant in the fc-th system. Moreover, g is the light-atom coupling strength and ^(^(tJo) the photonic (atomic) 
frequencies respectively (ft = 1 throughout the paper). The H^'^'"^ + H"'^* Hamiltonian can be diagonalized in a basis 
of mixed photonic and atomic excitations, called polaritons. On resonance between atom and cavity, the polaritons 
are created by operators P^ '" = | n±)j, {g,0\. The states | n±)j, = (| 5, n) j, ± \e,n— l)j.)/-\/2 are the polaritonic 
states (also known as dressed states) with energies Ei^ — nojd ± g\/n and | n)^, denotes the n-photon Fock state of the 
fc-th cavity. 

It has been shown that in an array of these atom-cavity systems the addition of a hopping photon term 
oc ^jictjaj+i + aja,_^-^), leads to a polaritonic Mott phase where a maximum of one excitation per site is allowed [5|]. 
This originates from the repulsion due to the photon blockade effect [2l[. In this Mott phase, the system's Hamiltonian 
in the interaction picture results 

Hi=jj:{pt'''p^^'+pt''p^t'')^ (3) 

k 

where J is the coupling due to photon hopping from cavity to cavity. Since double or more occupancy of the sites 
is prohibited; one can identify P^ '' with cr^ = f jj + icr^, where cr^. a^ and cr| stand for the usual Pauli operators. 
The system's Hamiltonian then becomes the standard XY model of interacting spin qubits with spin up/down 
corresponding to the presence/ absence of a polariton [5|. 

Let us now consider a linear chain of three coupled cavities with the two extremal ones doped with a two level 
system as shown in Fig. 1(a). Alternatively, as the central cavity in any case is undoped, one can simply replace it 
with an optical fiber of short length (so that the distance is greatly increased but the fiber still supports a single mode 
of frequency near those of the two cavities), which simplifies the setting even further, as shown in Fig. 1(b). For the 
purposes of description, we will use the three cavity setting remembering that everything applies to the case of two 
cavities linked by a fiber. The fact that a classical field can drive (i.e., pump energy into) the central cavity in a three 
cavity setting (as also shown in Fig. 1(a)) is replaced in the fiber setting by a coupler feeding light into the cavity (as 
also shown in Fig. 1(b)). 

Let tjj = |1— )j((7, 0| be the polaritonic spin operators for the end cavities (the index j = 1,2 labels the two end 
cavities) and a, a^ the field operators of the central empty cavity. Since the latter is not doped, there the field 
operators play the role of polariton operators and they couple to polariton operators of the ends cavities. Moreover, 
assuming that the central cavity (or fiber) is driven, the Hamiltonian describing the system dynamics will be 

2 
H = J^ (aja'^ + crjaj - Aa^a + aa'^ + a* a, (4) 

i=i 

where A = (^mid — ^poi is the detuning between the central cavity mode of frequency tOmid and the polaritons 
frequency uipoi = ujq — g. Furthermore, a is the product of the coupling of the driving field to the central cavity 
field (say G) and the amplitude of the driving radiation field (say a). We also assume that A is much smaller than 
the atom-light coupling in each of the outer cavities, so that only the ground level \g) = \g, 0) and first excited level 
|e) = {\g, 1) — |e, 0))/-\/2 of the polaritons are involved (i.e., the polaritons are still good as qubits). 

Suppose that the polaritons decay with the same rate 7 (this is the effective decay rate of the polariton due to both 
the decay of the cavity field and the atomic excited state), and the cavity radiation mode with rate k. The quantum 
Langevin equations describing the dynamics will be |26l | 



'3 



iJaa] - ^a, + ^of , j = 1 , 2 (5) 

a = iAa — iJ (a\ + a^) — ia ^ na + \/2kci™ , (6) 
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FIG. 1: The system under consideration, a) The cavities are coupled through direct photon hopping, b) The cavities are 
coupled through a fiber. The extremal cavities in each configuration are interacting with a two level system that could be an 
atom or a quantum dot depending the implementation technology used, c) The photon blockade allows for the ground and 
first dressed states of each atom-cavity system to be treated as a (polaritonic) qubit. 



where the superscript in denotes the vacuum noise operators. 

If K 3> J the radiation mode can be adiabatically cUminatcd in such a way that 
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(7) 



Moreover, if the quantities J /{2^/k) and a/{2^/K) are large compared to the amplitude standard deviation of the 
fluctuating vacuum field, the last term in Eq.([71) can be neglected and 



A- 



(ai + era) + 



A- 



Inserting Eq.® mto Eqs.®, we get[2£ 
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corresponding to an effective Hamiltonian for polaritons of the type 
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We are using ^ and 3 to denote the real and imaginary part respectively. 

The dynamics of the polaritons can now be described by the master equation [26 



2 ^ 



(8) 

(9) 
(10) 

(11) 



(12) 



where Lj — y^2r/aj are the Lindblad operators. 



III. STEADY STATE ENTANGLEMENT 



At the steady state Ea. (fT2)) becomes 



= -iC 



aial 



0-1^2, /9 



-^C 



al,p 



-iC Wl +cr2,p] 



-2aipal — (j\(7ip — po-jcri + 2a2p(72 — <^2'^2P — P'^2'^2, 



(13) 



where C = 3ff[J^/7(A + m)] and ^ == aJ/"f{A + in). 

The steady state sohition of Eq. (|13p can be found by writing the density operator and the other operators in a 
matrix form, in the basis B = {|e)i|e)2, |^)i|e)2, |e)i|.g)2, |.9)i|5)2}- Let us parametrize the density operator as 



P = 



A Bi+ iB2 Ci + iC2 Vi 

Bi-iB2 E Ti+iT2 Gi 

Ci — iC2 T\ — 1^2 H Xi 

Vi~iV2 gi-iV2 Ti-iT2 \-A 




(14) 



where the matrix elements also respect the requirement that Tr{p} = 1. The matrix representation of the other 
operators comes from 
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(15) 



By using matrices ((T4| and (fT5|) in the r.h.s. of Eq. p^ . we get a single complex matrix M which must be equal to 
zero. Then, equating to zero the entries of M we get a set of equation for the entries of p. Since M is Hermitian we 
can consider 



Mjj = 0, i,k = 1,2,3,4 



3?{M,fe} = 0, fc>j 
3{Af,fe} = 0, fc>j 



(16) 
(17) 
(18) 





-6-4 - 3Bi 


- CC2 + i^vx - ei2?2 + ii£ 




a^ - 302 


+ cci + iiVx + i2V2 - iiE 




-iiA-iB2- 


- 3Ci + 6^1 - iiV2 + 6^1 




iiA + C61 - 


- 3C2 + ixv^ + i2V2 - iiTr 


6^2- 


- ixB2 - aci 


-eiC2- 2151+6^1 +6^2 


6fii- 


- iiB2 + CiCi 


-i2C2-2V2-ixQx^ilQ2 



so to have a set of 16 linear equations. They are not all independent because of the 15 unknown parameters 
(A,B\^B2,C\,C2,T>\,'D2,£,T\,T2,Q\,Q2^'H^T\,X2). Explicitly the set of equations results 

- 4^ + 26Si - 2^162 + 26Ci - 2aC2 - 0, 

-6-^i-a-^2 = 0, 
a-^i-6-^2 = 0, 

viiT2^i2n = 0, 
K2^2-eiw = 0, 

f 62:1 + ^2:2 = 0, 
-6X1 + 62:2 = 0, 

2^ - 26^1 + 26^2 - 2£ - 2C.F2 + 26^1 - 26^2 = 0, 

-6^1 + 6^2 - 6C1 + 6C2 - 2.F1 + 6^1 - 6^2 + 62:1 - ^ii2 = 0, 

6S1 + 6^2 - 6C1 - 6C2 + C^ - 2T2 + 6^1 + 6^2 -cn- cai - 62:2 = 0, 

-6^ + 2Ci - 62?i + 6^2 - 26^ - 6-^1 - 6-^2 -Gi- ^iH + CI2 = -6, 

6^ + 2C2 - iiVi - 6i?2 + 26^ + 6-^1 - 6-^2 -52 + 6W - ^i = 6, 

2yt - 26Ci + 26C2 + 2CJ^2 - 2W + 26X1 - 26^2 = 0, 

-6^ + 2Bi - 6:Pi + 6^2 - 6^ - 6-^1 + ei-^2 + (02 - 26^ - ii = -6, 
6^ + 2B2 - iiVi - 6^2 + 6f + 6-^1 + 6-^2 - CQi + 26^ -^2 = 6, 

2^-26^1 + 26^2 + 27^-262:1 + 262:2 = 0, (19) 



where 6 = ^{i} and 6 = ^{C}- 




FIG. 2: Concurrence C versus (^2in(r§ 
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Solving analytically the above set of equations we obtain for ^2 
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where 



2\2 



C' + (1 + 2ef ) 



(20) 



(21) 



Notice that for <^i = we have formally analogous solutions that lead to the same physical result, hence they are not 
reported. 

Now that we know the stationary density matrix, we can use the concurrence as measure of the degree of entangle- 
ment [H 



C{p) = max {0, Ai — A2 — A3 — A4} , 
where A^'s are, in decreasing order, the nonnegative square roots of the moduli of the eigenvalues of pp with 

p^ (crfcr|)p*(crfcr|), 
and p* denotes the complex conjugate of p. With respect to the basis B it results 



P = 



1-A-S-H -I1-1I2 -Qi- 1Q2 'Di + 1V2 

-Xi + 1X2 'H Ti+ iT2 -Ci - iC2 

-Q\ + iQ2 J^i - iJ^2 £ -Bi - 162 

V1-W2 -Ci+iC2 -Bi+iB2 A 



(22) 
(23) 

(24) 



In Fig[2] we show the concurrence as a function of C and ^1 (the cases ^ = and ^ = give the same numerical 
results for the concurrence). Notice that by increasing (^, the concurrence increases quite slowly, and a maximum 
amount of entanglement is approximately 0.3 for C = 10 and ^1 = 2.135. This is similar to the amount of stationary 
entanglement achievable with an effective interaction of the kind erf erf when combined with an intricate feedback and 
cascading p7| . 

One could try to employ entanglement witnesses to detect this entanglement [25l |. A witness can be constructed 
from the density matrix corresponding to the maximum value of the concurrence. This would be a traceclass operator 
W in the Hilbert space of the two polaritonic qubits such that Tr[Typ] > for all separable states while Tr[VF/9] < 
for the considered entangled state. The form of such a witness in the Pauli decomposition results 



j,k—id,x,y,z 



(25) 



with a^'^ = I. In Fig[3] we show the coefficients Cj^k for the entanglement witness coming from the density matrix 
corresponding to the maximum value of the concurrence in Figl^l As we can see, the elements with the most 



FIG. 3: 
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significant weights (greater than 0.05) for measuring the witness, correspond to total of five measurements: two 
separate measurements of a^ in each polariton, and two joint measurements trf (E) erf and af (E) af- 

The values of C and ^i used in FigOto get maximal entanglement would correspond to A = 10 J, k = 10 J, G = 7 = 
O.OIJ and the pumping coherent field was also taken to have roughly a hundred photons. J is tunable and depends 
on the coupling of the photonic modes between neighboring cavities. Assuming this to be of the order of lO^'^Hz, this 
would correspond to a cavity dissipation rate k w lO^^Hz and a polaritonic decay rate 7 w lO^Hz. These correspond 
to 0.1 nanoseconds lifetime of the cavity field and to ten nanoseconds for the polaritonic excitations at the two ends, 
which are within the near future in technologies like coupled toroidal microcavities and coupled superconducting 
qubits [221 . Coupled defect cavities in photonic crystals arrays are also fast approaching this dissipation regime and 
are extremely suited in fabrication of regular arrays of many coupled defect cavities interacting with quantum dots 
[23|. In all technologies, an increase in J, in coupling between the cavity modes, the requirements on the various 
lifetimes of the polaritonic and photonic field modes can be further reduced. 



IV. CONCLUSION 



To summarize, this paper presents an example of entangling two qubits in the presence of dissipation despite the 
fact that each qubit has a continuously decaying state. The entanglement is not transient but stationary, and thereby 
easy to verify in an experiment, for which there is also a relevant witness. Though the amount of entanglement is 
not maximal, it is still very interesting as it is for a completely open system. As opposed to the typical case of, say, 
many-body systems or even the case of two purely atomic qubits in a single cavity or extremely close as to be able to 
directly interact, here there is the added advantage that the entangled qubits are easily individually accessible (being 
encoded in distinct atom-cavity systems) for measurements. It is worthwhile to point out an existing scheme to have 
steady state entanglement between entities in distinct cavities entangles atoms [l9[ (as opposed to polaritons) and is 
much more intricate. 

It is very interesting and counterintuitive that only a classical laser field driving the central cavity/connecting fiber 
was necessary to entangle the polaritonic qubits. A scheme feasible with current or near future technology and able to 
verify polaritonic entanglement as the one we have suggested in this paper, would be a significant first step towards the 
realization of the plethora schemes to simulate many-body systems and quantum computation using coupled cavities. 
Moreover, the model would also deserve to deepen counterintuitive properties of entanglement against noise (see e.g. 
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